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Abstract Heterogeneous A:-core percolation is an extension of a percolation model 
which has interesting applications to the resilience of networks under random dam- 
age. In this model, the notion of node robustness is local, instead of global as in 
uniform A:-core percolation. One of the advantages of A:-core percolation models is 
the validity of an analytical mathematical framework for a large class of network 
topologies. We study ternary mixtures of node types in random networks and show 
the presence of a new type of critical phenomenon. This scenario may have use- 
ful applications in the stability of large scale infrastructures and the description of 
glass-forming systems. 



1 Introduction 

Percolation, with its many modifications and extensions, is a problem with a ven- 
erable past and many applications in the most diverse disciplines [18]. With the 
rapid development of network science, several new problems and models have been 
introduced. Among the different directions, fe-core percolation constitutes a devel- 
opment which, in spite of its somewhat simple definition, is able to encapsulate a 
number of interesting problems which can often be approached with a robust math- 
ematical formalism. Given a network, a fe-core is defined as the subnetwork where 
each node has at least k neighbours in the same subnetwork. A £-core can equiva- 
lently be defined as the subnetwork remaining after a culling process consisting in 
recursively removing all the nodes with degree lower than k. fe-core percolation has 
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applications in many different disciplines including jamming [15], neural networks 
[8], granular gases [1], evolution [14], social sciences [13] and the metal-insulator 
transition [4]. The application of fe-core percolation we are interested in here, how- 
ever, refers to the stability of a network under random damage. A fraction ( 1 — p) of 
nodes is removed (together with the adjacent edges) and the size M# of the largest 
&-core cluster is studied as a function of p. The fe-core strength can vanish ei- 
ther continuously or discontinuously, a scenario that can be described in terms of 
phase transitions (second and first order, respectively) and can be important in the 
stability of large scale infrastructures [9, 10]. A recent extension of fc-core percola- 
tion, named heterogeneous &-core (HKC) percolation considers the threshold fc, as a 
local property, and allows therefore the mixing of different critical phenomena. An 
analytical formalism has been introduced to approach this model on networks and 
binary mixtures of thresholds k have been investigated [2, 5, 6]. On a wide class 
of networks, this model displays phase diagrams which are topologically equiva- 
lent to the ones calculated for a recently introduced spin model of glass-forming 
systems [17, 16]. This model can be seen as a heterogeneous development of the 
Fredrickson-Andersen (FA) model of facilitated spins, where the facilitation con- 
sists of a local constraint in the number of spins down in order for the considered 
spin to be able to flip [11, 17]. It has been shown that this model reproduces char- 
acteristic signatures of glass-forming systems called glass-transition singularities, 
which correspond to distinctive critical phenomena for the appropriate choice of 
the parameters [12, 16]. One of the glass-transition singularities which has not yet 
been explored is the so-called A4 singularity, characterised by the coalescence of 
two critical points. As the HKC model appears to reproduce all the relevant criti- 
cal phenomena for binary mixtures, it is interesting to look for this singularity as 
well. The purpose of this paper, therefore, is to investigate, for the first time, ternary 
mixtures in HKC percolation and calculate the critical point corresponding to an A4 
singularity. 

In Section 2 we present the analytical formalism used here (from [2]), in Sec- 
tion 3 we sketch the behaviour of the model for binary mixtures, in Section 4 we 
show a case of ternary mixture with an A4 singularity and in Section 5 we give the 
conclusions. 



2 Formalism 

Let k{ be the A:-core threshold of the node ;. A fraction (1 — p) of nodes are randomly 
removed: the problem consists in calculating the size of the HKC (if it exists). We 
consider the configuration model of random networks, defined as the maximally 
random network with a given degree distribution P(q). The configuration model has 
the property of being locally tree-like, i. e. the number of finite loops vanishes for 
infinite networks. This property allows to consider, in an infinite network, the HKC 
equivalent to the (fc, — l)-ary subtree, defined as the tree in which, as we traverse it, 
each encountered vertex has at least k t — 1 child edges. Then we can write a self- 
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consistent equation for Z, the probability that a randomly chosen node is the root of 
a (k { - l)-ary subtree [2]: 

*-"i^ l l( , 7>-^ + 



q=k b (?) l=k b -\ 

+p{\-r-s) £ ^ £ f'T^d-Z)^-'. (1) 

where the three thresholds k = (k a ,kb,k c ) are randomly assigned to nodes with prob- 
ability r, s and (1 — r — s), respectively. In this paper we always assume k{ > 2, so 
we do not need to consider the case where there may be finite clusters in the HKC. 
Due to the absence of finite loops, no finite HKC clusters can exist if ki > 2 for all 
nodes i [10, 2]. 

We can then write the probability M a j, c that a randomly chosen node is in the 
HKC, for a mixture of three types of nodes: 



M abc (p) =Pr£ P{q) t (f)z'(l -Z)"- l + 

q=k a l=k a V/ 
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3 Binary mixtures in heterogeneous fc-core percolation 

In homogeneous &-core percolation it is known that, for networks with a fast de- 
creasing degree distribution, namely P(q) < 1 /q r with y > 3 for q — > a randomly 
damaged network collapses continuously for k < 2 and discontinuously for k > 3 
[7, 3, 10]. In recent papers, it has also been shown that in the case of binary mix- 
tures HKC percolation is characterized by a few different critical phenomena. For 
example, in mixtures where the two values of k are associated to continuous and 
discontinuous transitions, respectively, a critical point (also called A3 singularity or 
cusp singularity) is usually observed [2, 6]. That is the case, for instance, of the 
mixtures k = (1,3) or k = (2,4), where the line of first order transitions ends in a 
critical point and a line of second order transition intersects the former. The case 
k = (2,3), though, is quite peculiar as the critical line exactly matches the line of 
first order transitions giving rise to a tricritical point (Fig. 1) [5]. 
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Fig. 1 Comparison of the phase diagrams of the cases k = (1,3), k = (2,3) and k = (2,4) for 
Erdos-Renyi graphs with Zi = 10. Continuous lines represent second order phase transitions, 
whereas dashed lines represent first order phase transitions. The red dot represents a tricritical 
point; the other dots are critical points. In terms of the notation given in Section 2, a binary mixture 
corresponds for example to s = 1 — r. 

When both values of k are characterised by a first order transition, instead, the 
phase diagram displays a critical point only when the two values of k are different 
enough [6]. Fig. 2, for instance, shows that the mixture k = (3, 8) has a critical point 
whereas k = (3,4) does not have one. 



4 The A4 singularity 

From the analogy between HKC percolation and facilitated spin models, it is in- 
teresting to investigate whether the singularities observed in models of the glass- 
forming systems are also present in HKC percolation. In particular, the fact that all 
the critical phenomena observed so far in HKC percolation are in the same class 
of universality as the ones observed in the FA model with heterogeneous facilita- 
tion is quite remarkable [16, 6]. The MCT of glass-forming systems also predicts 
an "A4 singularity", meaning the coalescence of a critical point into a line of first 
order phase transitions [12]. The A4 singularity is also named a swallow-tail bifur- 
cation. As this scenario is associated with a three-parameter theory, it is reasonable 
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Fig. 2 Phase diagram of the cases k = (3, 8) and k = (3,4) for the Erdos-Renyi graph (zi = 30). 
Symbols are as in Fig. 1. While for k = (3,4) there is only one line, in the case k = (3,8) there are 
two lines of discontinuous transitions and a critical point. 



to expect that the corresponding singularity in our HKC model can be observed in 
a three-component mixture. Moreover, we have seen in Section 3 that binary mix- 
tures with ki > 3 are either characterised by a single line of first order transitions or 
two lines of first order transitions with a single critical point. It is therefore natural 
to consider a ternary mixture which interpolates between the two regimes. We now 
show that an A4 singularity is indeed present in the ternary mixture k = (3,5, 8) on 
Erdos-Renyi graphs. 

Equation (1) can be re- written as 



/358(Z) = p 



(3) 



where 



rziZ+s £(^ +(1 



11=1 



n 



n 



(4) 



At every fixed fraction of remaining nodes p, depending on the values of the param- 
eters r and s, /35s has one or two maxima in Z. If the second maximum is higher 
than the first one, there is only one first order transition between a percolating and 
a non-percolating HKC. If the first maximum is higher than the second one, there 
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are two first order transitions: one between a high-fc and a low-fc phase and another 
towards a collapsing HKC. Typically, the second maximum disappears at a critical 
point, which can be found by imposing the condition: 

/^ 58 (Z)=/^ 8 (Z)=0. (5) 

This defines a locus in the (r, s) plane corresponding to a line of critical points. At 
high values of s, however, the critical point disappears and only a single line of 
first order transitions survives. This is due to the fact that the 5-nodes interpolate 
between the other two values of k and there is no critical point in either the mixture 
k = (3,5) or k = (5,8). The condition of an A4 singularity corresponds to the onset 
of this behaviour and is defined by the condition: 

f 35S (Z)=fl 5S (Z)=f£ s (Z)=0, (6) 

which yields 

r* =0.3140687639806 s* = 0.1831697392197 (7) 

which is the position of the A4 singularity on the phase diagram. Fig. 3 shows the 
phase diagram in the plane (p, r) for a few values of s. 

The critical exponent j3 is defined by the vanishing of the order parameter in 
approaching a phase transition: M^%{p) — M^ss(Pc) ~ (p — Pc)^ ■ As in binary 
mixtures, the critical exponent at the discontinuous phase transition is J3 = 1 /2, 
whereas it becomes j3 = 1/3 at the critical point. In the vicinity of the A4 singular- 
ity, we have (p — /?* ) <~ /35s (Z* ) — /35s (Z) <~ (Z — Z* ) 4 , from which it follows that 
■^358 (p) — ^358 (p*) ~ (j> — P*) 1 / 4 ) asM35g is linear in p and Z in the vicinity of the 
transition. Thus, the values of the critical exponent j3 can be summarized as follows: 

{1/2 hybrid transition r/r c ,i / 
1/3 critical point r = r Cl s^ s* (8) 
1 /4 A4 singularity r = r* , s = s« 

The unique value of /3 = 1/4 identifies the critical point as a A4 singularity. This phe- 
nomenon has an interesting interpretation. In a networks formed by types of nodes 
of very different fragilities (as in k a = 3, k c = 8), the collapse due to random damage 
can be catastrophic, with the presence of multiple discontinuous transitions. The in- 
troduction of nodes of intermediate fragility (kf, = 5) can weaken the critical point up 
to disappearance, making a significant region of the phase diagram stable (Fig. 3). 
This purely mathematical remark could be useful in principle in planning the stabil- 
ity of large scale infrastructures. Another interesting application is in modelling the 
glass transition. The presence of A4 singularities is predicted by the mode-coupling 
theory of glass-forming systems [12]. Due to the strong similarities between HKC 
percolation and facilitated spin models, it is reasonable to expect that this singular- 
ity should also be present in suitable ternary mixtures of facilitation parameters in 
such spin models [17, 16]. 
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Fig. 3 Phase diagram of the mixture k = (3,5,8) for Erdos-Renyi graphs (zi = 30). r and s are 
the fractions of nodes of type 3 and 5, respectively, and p c is the fraction of undamaged nodes at 
a phase transition. Lines of discontinuous transitions are plotted in the plane (p, r) at three fixed 
values of s. In the absence of type 5 (s = 0), there are two lines of first order transitions. One of 
them ends in a critical point (black dot). At the critical fraction s», the two lines touch at a A4 
singularity (red star). For s > s*, there is only a single line of discontinuous transitions. 



5 Conclusions 

In this paper we have applied the formalism of heterogeneous £-core percolation to 
ternary mixtures of thresholds k. In particular, we have calculated the phase diagram 
of a ternary mixture for Erdos-Renyi graphs. This mixture displays a characteristic 
A4 singularity, i. e. a critical phenomenon characterised by the merging of a critical 
point with a distinct line of discontinuous transitions. A peculiarity of this point 
is the change in the critical exponent j3 of the order parameter, which uniquely 
assumes the value 1/4. The characteristics of the studied phase diagrams may give 
useful information in designing large scale infrastructures to be resilient to random 
damage. This singularity is also predicted by models of the glass transition [16] and 
our model appears to be in the same universality class of kinetic spin models with 
heterogeneous facilitation. 
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